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ABSTRACT 

The whole sky differential star counts (DSC) with 1 degree resolution are retrieved 
from 2MASS online data service. Galaxy with double exponential thin and thick disks 
and a single power law luminosity function (LF) is used to interpret the 2MASS data. 
The slope of the DSC appears roughly isotropic over the whole sky, the average value 
is ~ 0.32, which corresponds to a power law index ~ 1.8 of the LF. We find that the 
scale-length and scale-height the thin disk are ~ 3.0 kpc and ~ 245 pc, and those of 
the thick disk are ~ 3.0 kpc and ~ 780 pc. The ratio of the thick disk to the thin disk 
is ~ 7%. The location of Sun above the disk is ~ 15 pc. A comparison of the data and 
model and their discrepancy are also provided. 
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1 INTRODUCTION 

Star counts provides us a good way to investigate the struc- 
ture of the Milky Way. In 1980 Bahcall & Soneira proposed 
a two component model, disk and halo, to interpret Galac- 
tic structure (Bahcall & Soneira 1980). However, in 1983 
Gilmore & Reid identified the disk could be further divided 
into two parts, thin and thick, for better fitting (Gilmore 
& Reid 1983). Since then, many works based on different 
observations with limited or specific sky area have been car- 
ried out (Karaali et al. 2004) and indicate that the thick 
disk is an individual component. The scale-heights of the 
thin disk in these works are similar, ranging from 240 kpc 
to 350 kpc. However, the scale-height of the thick disk have 
a much larger range and tends to have smaller value with 
larger local stellar number density (Chen et al. 2001). The 
Two Micron All Survey 2MASS (Skrutskie et al. 2006), pro- 
vides an opportunity to use large sky coverage star counts 
to investigate the Galactic structure and more detail struc- 
tures were recognized. The warp found on 2MASS is the 
same as that from gas and dust. Also, the scale- height of 
the thick disk increases with the Galactocentric radius for 
R > b kpc (Cabrera-Lavers et al. 2007). However, the best 
model to describe the major smooth components, thin and 
thick disks, of the Milky Way has not been fixed yet. 

In this contribution, we use the whole sky 2MASS point 
source catalog (hereafter, 2MASS PSC) to derive the overall 
structural parameters and the luminosity function of the 
Milky Way. 
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Figure 1. Left:The Ks band DSC of 2MASS (black plus) and 
model (blue cross) at [l,b] = [270,10]. Right After removing the 
slope, DSC clearly shows a dip. The black plus is 2MASS data, 
and the blue cross is model. 

2 2MASS POINT SOURCE CATALOG 

Sources with SNR ^ 5 and detected in all J, H, Ks bands 
were selected from 2MASS PSC. Its completeness rate cor- 
responding lOcr SNR limiting magnitude is 99%. By Hierar- 
chical Triangular Mesh, the whole sky is divided into 32768 
nodes with ^ 1 degree resolution. The data was retrieved 
directly from 2MASS online data service via virtual obser- 
vatory protocol. We only use Ks band data in this work. Due 
to the small number of stars brighter than 8 magnitude and 
incompleteness fainter than limiting magnitude, we analyze 
data between 8 to 14 Ks magnitude only. The differential 
star counts (hereafter DSC) is well approximated by a sin- 
gle power law in this region (see Fig.l left). 



3 THE MODEL AND ITS FEATURES 

Star count depends on the distribution of stars in space and 
luminosity, i.e., stellar number density profile of the Galaxy 
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All Sky 2MASS Slope 




Figure 2. The whole sky 2MASS DSC slopes in magnitude 8 to 
14 seem homogeneous. The average value is ~ 0.32. 




Figure 3. The method of maximum likelihood is used to find hr 
and hz. The hr and hz of the thin disk are ~ 3.0 kpc, ~ 245 pc 
and that of thick disk are ~ 3.0 kpc, ~ 780 pc. The ratio of the 
thick to the thin disk is ~ 7%. 

and local luminosity function. The density profile, n(/ir, hz)^ 
in our model contains a thin disk and a thick disk. Both are 
exponentially decayed along the Galactocentric radius, i?, 
and the distance above or below the Galactic plane, Z, 

n{r, z) = no exp ^ j— j exp ^ j-— j , (1) 

where hr and hz are the scale-length and scale-height of the 
disk, no is the local density in the solar neighborhood, Rq 
and Zq are the sun location. In our model the luminosity 
function (hereafter LP), < L < L2) oc is a single 

power law, where Li, L2 are lower and upper cutoffs. More- 
over, we assume that the LP does not depend on location. 
The synthetic DSC can be divided into 3 regions, bright, 
modest and faint, with slopes 0.6, 0.4(7 — 1) and a cutoff, 
respectively (See Pig.l left). Therefore, we can derive 7 from 
the slope in the modest region. Moreover, the turning point 
and cut off can be determined by hr, hz, Li and L2. The 
New COBR/IRAS extinction model (hereafter NCI model) 
was adapted for Galactic latitude \b\ ^ 10° area (Chen et 
al. 1999). 




Figure 4. The number ratio of the 2MASS data to the model. 
Left: towards anti-Galactic center. Right: towards the Galactic 
center. 



4 RESULTS AND DISCUSSION 

4.1 Luminosity Function 

The differential star counts between Ks magnitude 8 to 14 
was used to derive the slope of the DSC of the whole sky. 
The slope seems roughly isotropic over the whole sky, and 
the average value is ~ 0.32, and hence the power law index of 
the LP, 7 ~ 1.8 (see Pig. 2). If the slope is removed from the 
DSC, a small but significant dip appears (see Pig.l right). 
Although the depth and width of the dip are hard to quan- 
tify, it can only be reproduced by some deficiency on the LP. 
We strongly believed that the dip is not caused by the fluc- 
tuation on density profile. In this case, the deficiency on LP 
can help us to quantify hr and hz , and vice versa. According 
to the derived hr and hz of the thin disk, the deficiency on 
the LP is located at Ks absolute magnitude. 

4.2 Density Profile 

The method of maximum likelihood was applied to find the 
best fit hr and hz- In region with \b\ ^ 30°, the thin disk hr 
and hz are ^ 3.0 kpc, ^ 245 pc, and those of the thick disk 
are ~ 3.0 kpc, ~ 780 pc. The ratio of the thin to the thick 
disk is ^ 7% (see Pig. 3). Zq is derived from the number 
ratio of north to south Galaxy, and the value is ^ 15 pc. 

4.3 Comparison Between Data and Model 

Pig. 4 shows the number ratio of the 2MASS data to the 
model. Since we do not include the Galactic bulge in our 
model, there is an obvious discrepancy at the Galactic cen- 
ter. Our model generates less star in the high latitude area. 
This discrepancy can not be explained away by adding thick 
disk in our model. In addition, there is another discrepancy 
in the anti-GC direction. It might be a result of a warp or 
an increasing hz with Galactocentric radius (Cabrera-Lavers 
et al. 2007). The NCI model has large uncertainty in low 
Galactic latitude area, therefore our model produces unrea- 
sonably large star number in these areas. If we believe that 
stellar distribution on the disk is smooth, the star count can 
constrain the extinction model. 
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